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ABSTRACT 

We present the first weak gravitational lensing analysis of the completed Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS). We study the 64 deg 2 WI field, the largest of the CFHTLS- 
Wide survey fields, and present the largest contiguous weak lensing convergence "mass map" yet 
made. 

2.66 million galaxy shapes are measured, using a Kaiser, Squires and Broadhurst (KSB) pipeline 
verified against high-resolution Hubble Space Telescope imaging that covers part of the CFHTLS. 
Our i'-band measurements are also consistent with an analysis of independent r'-band imaging. The 
reconstructed lensing convergence map contains 301 peaks with signal-to-noise ratio v > 3.5, consistent 
with predictions of a ACDM model. Of these peaks, 126 lie within 3'.0 of a brightest central galaxy 
identified from multicolor optical imaging in an independent, red sequence survey. We also identify 
seven counterparts for massive clusters previously seen in X-ray emission within 6 deg 2 XMM-LSS 
survey. 

With photometric redshift estimates for the source galaxies, we use a tomographic lensing method 
to fit the redshift and mass of each convergence peak. Matching these to the optical observations, we 
confirm 85 groups/clusters with Xreducod < 3.0, at a mean redshift (z c ) = 0.36 and velocity dispersion 
(a c ) = 658.8 km s" 1 . Future surveys, such as DES, LSST, KDUST and EUCLID, will be able to 
apply these techniques to map clusters in much larger volumes and thus tightly constrain cosmological 
models. 

Subject headings: cosmology: observations, - galaxies: clusters: general, -gravitational lensing: weak, 
-X-rays: galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies are the largest gravitationally 
bound structures in the universe. The number and mass 
of the biggest clusters are highly sensitive to cosmological 
parameters including the mass density fl mi the normal- 
ization of the mass power spectrum erg (e.g., Press & 
Schechter 1974; Frenk et al. 1990; Eke et al. 1996; Sheth 

1 Department of Physics and Tsinghua Center for Astro- 
physics, Tsinghua University, Beijing, 100084, China 

2 Laboratoire d'Astrophysique de Marseille, CNRS-Universite 
de Provence, 38 rue Frederic Joliot-Curie, 13 388 Marseille 
Ccdex 13, France 

3 Centre de Physique des Particules de Marseille, 
CNRS/IN2P3-Luminy and Universite de la Mediterranee, 
Case 907, F-13288 Marseille Cedex 9, France 

4 Department of Astronomy, Peking University, Beijing, 
100871, China 

5 Dark Cosmology Centre, Niels Bohr Institute, University 
of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, 
Denmark 

6 Institute for Astronomy, Royal Observatory, Blackford Hill, 
Edinburgh EH9 3HJ, UK 

7 California Institute of Technology, MC 350-17, 1200 East 
California Boulevard, Pasadena, CA 91125, USA 

8 Jet Propulsion Laboratory, California Institute of Technol- 
ogy, Pasadena, CA 91109, USA 

^ Canada France Hawaii Telescope, 65-1238 Mamalahoa Hwy, 
Kamuela, HI 96743, USA 

10 Department of Physics & Astronomy, University of Victo- 
ria, Victoria, BC, V8P 1A1, Canada 

11 National Research Council of Canada, Herzberg Institute 
of Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 
2E7, Canada 

12 Institude d'Astrophysique de Paris, UMR 7095, 98 bis 
Boulevard Arago, 75014 Paris, France 

' Email address: shanhuanySgmail . com 



& Tormen 1999), and the dynamics of dark energy (e.g., 
Bartelmann et al. 2006; Francis et al. 2009; Grossi & 
Springel 2009). Understanding the properties of clusters 
is vital to test theories of structure formation and to map 
the distribution of cosmic matter on scales of ~l-10Mpc. 

Theoretical predictions of structure formation deal di- 
rectly with the total mass of clusters; measurements 
are restricted to indirect proxies that can be observed. 
Contaminating the translation between theory and ob- 
servation are large uncertainties in the interpretation 
of galaxy richness, X-ray luminosity /temperature and 
the Sunyaev-Zeldovich decrement (e.g. Bode et al. 2007; 
Leauthaud et al. 2010). Weak gravitational lensing, the 
coherent distortion of galaxies behind a cluster, can po- 
tentially provide direct measurements of the total mass 
regardless of its baryon content, dynamical state, and 
star formation history. 

By measuring the shear (coherent elongation) of 
many background galaxies, we can reconstruct the two- 
dimensional (2D) weak lensing convergence map, which 
is proportional to density projected along each line of 
sight. Peaks in the convergence map with high signal- 
to-noise ratio v generally correspond to massive clusters 
(Hamana et al. 2004; Haiman et al. 2004). Since the 
three-dimensional (3D) shear signal should increase be- 
hind those clusters in a predictable way that depends 
upon only the lens-source geometry, we can also use pho- 
tometric redshift estimates of the background galaxies 
(from multi-band imaging) to measure the redshift and 
mass of each foreground cluster (Wittman et al. 2001, 
2003; Hennawi & Spergel 2005; Gavazzi & Soucail 2007). 
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Systematic weak lensing cluster searches have only re- 
cently become practicable. Miyazaki et al. (2002) used 
Subaru/Suprime-Cam in excellent seeing conditions to 
find an excess of 4.9 ± 2.3 convergence peaks with v > 5 
in an area of 2.1 deg 2 . Dahle et al. (2003) and Schirmer 
et al. (2003) each identified several shear-selected clus- 
ters with redshifts z ~ 0.5 determined from two-color 
photometry. Hetterscheidt et al. (2005) reported the de- 
tection of five cluster candidates over a set of 50 discon- 
nected Very Large Telescope/FORS images covering an 
effective area of 0.64 deg 2 , while Wittman et al. (2006) 
found eight detections in the first 8.6 deg 2 of the Blanco 
Deep Lens Survey. Gavazzi & Soucail (2007) presented 
a weak lensing analysis of initial Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS) Deep data covering 
4 deg 2 . They demonstrated that the image quality at 
CFHT is easily sufficient for cluster finding. Miyazaki 
et al. (2007) presented the first large sample of weak 
lcnsing-selected clusters in the Subaru weak lensing sur- 
vey, with 100 significant convergence peaks in a 16.7 deg 2 
effective survey area. Hamana et al. (2009) reported re- 
sults from a multi-object spectroscopic campaign to tar- 
get 36 of these cluster candidates, of which 28 were con- 
firmed (and 6 were projections along a line of sight of 
multiple, small groups). 

The main astrophysical systematic effect afflicting 
weak lensing cluster surveys is the projection of large- 
scale structure along the line of sight. Random noise 
is also added due to the finite density of resolved source 
galaxies and the scatter of their intrinsic shapes. Numeri- 
cal studies (White et al. 2002; Hamana et al. 2004) show 
that these contaminants significantly reduce the purity 
of cluster detection. To improve our analysis, we shall 
combine our weak lensing results with multi- wavelength 
imaging. Simultaneous detection of a weak lensing sig- 
nature plus an overdensity of galaxies with a single red 
sequence provides an unambiguous cluster identification. 
Furthermore, 3D lensing tomography using photometric 
redshifts from the multi- wavelength data can remove the 
other potential hurdles of: lensing signal dilution by clus- 
ter member galaxies, and identifying the redshift of weak 
lensing peaks when no corresponding galaxy overdensity 
is apparent. 

Here we present a weak gravitational lensing analysis 
of the 64 deg 2 CFHTLS- Wide Wl field, which is suffi- 
ciently large to contain several hundred galaxy clusters. 
Compared to the analysis of the CFHTLS-Deep survey 
by Gavazzi & Soucail (2007), our shallower CFHTLS- 
Wide imaging (and lower source galaxy density) will fa- 
vor the detection of higher mass, nearby clusters. The 
huge increase in survey area over any previous survey is 
expected to yield many more systems overall. In this 
paper, we shall primarily study the properties of the de- 
tected clusters, rather than the cosmology in which they 
are embedded. For this purpose, we adopt a default cos- 
mological model with Q m = 0.27, fi A = 0.73, cr 8 = 0.809, 
H a = 100 h km s" 1 Mpc~\ and h = 0.71. 

This paper is organized in the following way. In Sec- 
tion 2, we describe the CFHT and Hubble Space Telescope 
(HST) data used. In Section 3, we present the measure- 
ment of galaxy shapes in the CFHT imaging, and their 
calibration against measurements of the same galaxies 
in the HST imaging. In Section 4, we reconstruct the 



2D lensing convergence "mass map" signal, and extract 
a catalog of local maxima that represent cluster candi- 
dates. In Section 5, we search for optical counterparts of 
these candidates, dramatically cleaning the catalog. In 
Section 6, we investigate the full 3D lensing signal around 
each cluster, further cleaning the catalog when the lens- 
ing signal behind spurious peaks does not increase as 
expected with redshift - but obtaining an independent 
estimate of the cluster redshift when it does. We finally 
explore global scaling relations between cluster mass ob- 
servables, then conclude in Section 7. 

2. DATA 

2.1. CFHTLS-Wide T0006 imaging 

The CFHT Legacy Survey is a joint Canadian- 
French program to make efficient use of the CFHT 
wide field imager MegaPrime, simultaneously address- 
ing several fundamental questions in astronomy. Each 
MegaPrime/MegaCam image consists of an array of 9 x 4 
e2v CCDs with a pixel scale of 0".187 and a total field 
of view of ~ 1 deg 2 . The survey used most of the tele- 
scope dark and gray time from 2003 to 2008. We analyze 
CFHTLS-Wide imaging from the Terapix T0006 process- 
ing run, which is the first to include the complete survey 
and was publicly released on 2010 November 15 (Gora- 
nova et al. 2009). These data cover ~ 171 deg 2 in four 
fields (Wl, W2, W3, and W4) of which the 72-pointing, 
~ 64 deg 2 Wl field is the largest, and in five passbands 
(u', g', r' , i',and z 1 ) down to i ~ 24.5 and r ~ 25.0. 

Fu et al. (2008) showed that the z'-band exposures, 
taken in sub-arcsecond seeing conditions, provide the 
best image quality and resolve the galaxy population 
with highest median redshift. Resolving the shapes of 
more distant galaxies is vital for weak lensing analysis, 
since the strength of the shear signal is proportional to 
the ratio of the Lens-Source and Observer-Source dis- 
tances. We therefore choose to analyze the i'-band im- 
ages (mean seeing 0".73) in the contiguous Wl field. We 
also analyze the independent r'-band imaging to check 
the calibration of our shear measurements. We also use 
photometric redshift estimates for source galaxies ob- 
tained from the multicolor imaging (Ilbcrt et al. 2006; 
Coupon et al. 2009; Arnouts et al. 2010). 

Early releases of smaller regions of the CFHTLS have 
also been used to measure the weak lensing cosmic shear 
signal (Semboloni et al. 2006; Hoekstra et al. 2006; Fu et 
al. 2008). As the survey size has increased, the statistical 
errors have shrunk, and difficulty measuring shapes at 
a precision better than the statistical error has so far 
prevented publication of a cosmic shear analysis of the 
complete survey. However, weak lensing cluster searches 
are restricted by construction to regions of the survey 
where the signal is strongest, and the circular symmetry 
of our analysis removes the negative impact of additive 
shear measurement errors (cf. Mandelbaum et al. 2005). 

2.2. HST COSMOS Imaging 

The HST COSMOS survey (Scoville et al. 2007) is the 
largest contiguous optical imaging survey ever conducted 
from space. High resolution (0".12) imaging in the 
Ifsuw band was obtained during 2003-2005 across an 
area of 1.64 deg 2 that also corresponds to the CFHTLS 
D2 deep field. Any galaxies resolved by CFHT are very 
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Table 1 

SExtractor Configuration Parameters 



Pa,rcLmctcr 


Value 


DETECT.MINAREA 


3 


DETECT.THRESH 


1.0 


DEBLEND.NTHRESH 


32 


DEBLEND.MINCONT 


0.002 


CLEAN.PARAM 


1.0 


BACK.SIZE 


512 


BACK.FILTERSIZE 


9 


BACKPHOTO.TYPE 


local 


BACKPHOTO.THICK 


30 




36.800 36.600 36.400 36.200 36.000 



Righl ascension 

Figure 1. "Wl+2+3" pointing from the CFHTLS-Wide Wl field 
in i' band, showing masked regions. This pointing is representative 
of those with fairly poor image quality: the seeing of 0".78 is worse 
in only 24 of 72 (1 in 3) pointings. In our automated algorithm for 
masking diffraction spikes around bright stars, the basic shape of 
the star mask is predefined, and its size is scaled with the observed 
major axis of each star. 

easily resolved by HST, which therefore provides highly 
accurate shape measurements almost without the need 
for point-spread function (PSF) correction. We shall cal- 
ibrate our CFHTLS shape measurements against those 
from COSMOS by Leauthaud et al. (2010). 

Note that measurements of the shapes of individual 
galaxies from ground-based and space-based observations 
need not necessarily match exactly, even without shape 
measurement errors, because the different noise proper- 
ties of the data sets may make them most sensitive to 
different isophotes, which can be twisted relative to each 
other. The slightly different passbands may also empha- 
size different regions of a galaxy's morphology. However, 
across a large population of galaxies, these differences 
should average out, and a comparison of successful shear 
measurements between the two data sets should agree. 

3. GALAXY SHAPE MEASUREMENT 
3.1. Object Detection and Masking 



t — : — : t — : 1 1 r 




j i L 



18 20 22 24 26 

Mag 

Figure 2. Star selection (red points) in the planes of magnitude 
vs. flux radius (top) and magnitude vs. peak surface brightness 
(bottom). We find the latter more robust. The red points denote 
objects selected as stars for PSF modeling. The blue objects are 
spurious detection. 

We conduct shape measurement in both CFHTLS i' 
and r' bands. We detect astronomical sources in the im- 
ages using SExtractor (Bertin & Arnouts 1996). Our 
choice of the main SExtractor parameters is listed in Ta- 
ble [U and the data arc filtered prior to detection by a 3 
pixel Gaussian kernel. 

Near saturated stars, many spurious objects are found 
due to detector effects and optical ghosting. It would also 
be difficult to measure the shapes of real stars or galaxies 
in these regions, because of the steep background gradi- 
ents. We have developed an automatic pipeline to define 
polygonal-shaped masks around saturated stars, and all 
objects inside the masks are removed from our catalog. 
The masks in all images are then visually inspected; our 
automated pipeline fails in a few cases (mainly very satu- 
rated stars for which the centroid of the star measured by 
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Figure 3. Spatial variation of measured stellar ellipticities in the 
representative CFHTLS-Wide Wl+2+3 field, before (top) and af- 
ter (bottom) PSF anisotropy correction. The longest tick marks 
represent ellipticities of ~ 11%. The mean absolute ellipticity after 
correction is 0.62%. 

SExtractor was widely offset from the diffraction spikes) 
and those stellar masks are corrected by hand. An ex- 
ample of the masks for one CFHT pointing is shown in 
Figure Q] This pointing has slightly worse than average 
image quality, so we shall use it throughout this paper as 
a conservative representation of our analysis. After ap- 
plying all of our masks across the entire survey, the final 
effective sky coverage drops from 64 deg 2 to ~ 51.3 deg 2 
and 55.0 deg 2 for i' and r' bands, respectively. 

We shall employ the popular KSB method for galaxy 
shear measurement (Kaiser et al. 1995; Luppino & Kaiser 
1997; Hoekstra et al. 1998). In this method, the observed 
galaxy shape is modeled as a convolution of the (sheared) 
galaxy with the PSF, which is modeled as an isotropic, 
circular profile convolved with a small anisotropy. 

3.2. PSF Modelling 

To measure the shapes of galaxies, it is first necessary 
to correct them for convolution with the PSF imposed by 
the telescope optics and Earth's atmosphere. The chang- 
ing size and shape of the PSF across the field of view and 
between exposures can be traced from stars, which are 
intrinsically point sources. We identify stars from their 



e1 

Figure 4. Projection of the stellar ellipticities in the (ei , ez) plane 
before (black) and after (red) PSF anisotropy correction. The post- 
correction residuals are consistent with featureless white noise. 



constrained locus within the size-magnitude plane (Fig- 
ure (21(a)). Heymans et al. (2006) suggest using the full 
width at half maximum (FWHM). However, we find that 
FWHM is not robustly measured by SExtractor, so we 
instead use the /i max -magnitude plane (Bardeau et al. 
2005, 2007; Leauthaud et al. 2007), where /x max is the 
peak surface brightness (Figure [^b)). The red points in 
Figure [2] indicate the selected stars; our chosen locus re- 
flects a careful balance between obtaining sufficient stars 
to model the small-scale variations that we observe in the 
PSF pattern, and introducing spurious noise by includ- 
ing faint stars. The blue points are spurious detections 
of noise, cosmic rays, etc. (cf. Leauthaud et al. 2007). 

We then measure the Gaussian-weighted shape mo- 
ments of the stars, and construct their ellipticity. In 
addition to cuts in /x max and magnitude, we also exclude 
noisy outliers with signal-to-noise v < 100 or absolute el- 
lipticity e* more than 2a away from the mean local value, 
and we iteratively remove objects very different from 
neighboring stars. In 15 pointings with the worst im- 
age quality, including Wl+2+3, the PSF becomes larger 
than r g ~ 0".5 in the corners of the field of view, so we 
finally add these regions to the survey mask (and exclude 
galaxies in them from our weak lensing analysis). 

Having obtained our clean sample of stars, we con- 
struct a spatially varying model of the PSF across the 
field of view. In most pointings, we fit the ~ 30 stars in 
each of the 36 individual CCDs composing the MegaCam 
focal plane, using a polynomial of second order in x and 
y. For stacked data with large dithers, we use a higher 
order polynomial. Figures [3] and 0] show the stellar ellip- 
ticity before and after correction for the Wl+2+3 point- 
ing, using a weight function of default size r g to measure 
the PSF shape moments. The residual stellar ellipticity 
after correction is a consistent random scatter around 
zero, of width a ei ^0.01. 

The ellipticity of the PSF changes from the core to the 
wings. We measure the PSF shape using differently-sized 
weight functions and, when correcting each galaxy, use 
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Figure 5. Ellipticity of the PSF changes from the core to the 
wings. This shows the mean PSF ellipticity in the i'-band of the 
CFHTLS-Wide Wl pointing Wl+2+3 as a function of the size of 
the Gaussian weight function with which it is measured, before 
(black) and after (red) PSF anisotropy correction. The error bars 
show the rms scatter throughout that pointing. 

the same size weight function to measure both the PSF 
and galaxy shapes. Figure [5] shows the variation of mean 
stellar ellipticity as a function of the weight function size 
r g , before and after PSF anisotropy correction. 

3.3. Galaxy Shape Measurement 

Galaxies are selected as those objects with half light 
radius l.lr ; ^ SF < r/, < 4 pixels, where r^ SF is the size 
of the largest star, signal-to-noise v > 10, magnitude 
21.5 < i' < 24.5 and SExtractor flag FLAGS = 0. To also 
exclude blended or close pairs that could bias ellipticity 
measurements, we also cut objects with corrected elliptic- 
ity |e cor | > 1 and pairs of galaxies within 3". After survey 
masking and catalog cuts, the galaxy number density is 
n g ~ 11.5 arcmin -2 in an area of ~ 51.3 deg 2 of %'- 
band imaging; and n g ~ 7.9 arcmin~ 2 in A r i ~ 55.0 deg 2 
of r'-band data. Note that both are lower than the galaxy 
density n g ~ 38 arcmin -2 obtained in CFHTLS-Deep 
imaging by Gavazzi & Soucail (2007). This will restrict 
our detections to generally more massive clusters. Fig- 
ure [6] shows the magnitude distribution of the galaxies, 
and the redshift distribution of the 72 % (76 %) of galax- 
ies selected in the i' (r 1 ) bands that also have photometric 
redshift estimates by Arnouts et al. (2010). 

We then measure the shapes of all the selected galaxies. 
Our implementation of KSB is based on the KSBf9CE3 
pipeline (Heymans et al. 2006). This has been gener- 

14 http://www.roe . ac . uk/~heymans/KSBf 90/Home .html 
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Figure 6. Galaxy magnitude and redshift distributions. Top: 
number counts of galaxies per magnitude bin immediately after 
star-galaxy separation (solid histograms) and after all the lensing 
cuts (dotted histograms) for imaging in i' (black lines) and r' (red 
lines) bands. Bottom: the redshift distribution of galaxies used 
from the i' (black lines) and r' (red lines) bands. 

ically tested on simulated images containing a known 
shear signal as part of the Shear Testing Programme 
(STEP; Heymans et al. 2006; Massey et al. 2007) and 
the Gravitational lensing Accuracy Testing (GREAT08; 
Bridle et al. 2010) challenge. In all cases, the method was 
found to have small and repeatable systematic errors. 

If the PSF anisotropy is small, the shear 7 can be re- 
covered to first-order from the observed ellipticity e obs of 
each galaxies via 

/ psm \ 

7 = Pj 1 (> S - pi^* ) ■ (1) 

where asterisks indicate quantities that should be mea- 
sured from the PSF model interpolated to the position 
of the galaxy, P sm is the smear polarizability, and P 7 is 
the correction to the shear polarizability that includes 
smearing with the isotropic component of the PSF. The 
ellipticities are constructed from a combination of each 
object's weighted quadrupole moments, and the other 
quantities involve higher order shape moments. All def- 
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Figure 7. Choices for the way shear polarizability P 7 can be 
fitted to a galaxy population in the CFHT KSBf90 pipeline, to re- 
duce noise and bias in individual measurements. Linearly spaced 
contours compare our shear measurements of galaxies in a sub- 
set of the CFHTLS-Deep imaging, stacked to the depth of the 
CFHTLS-Wide survey, against measurements from the Hubble 
Space Telescope. Dashed lines show the best-fit relation -yCFHT = 
(m— l)7 1 iST +c. The four panels illustrate various fitting schemes. 
Top-left: raw (noisy) P 7 measurements from each galaxy, without 
fitting. Top-right: fitted as a polynomial in galaxy size P 7 (r^). 
Bottom-left: fitting function P 7 (rj,, mag) from Fu et al. (2008). 
Bottom-right: best-fit rational function P 7 (r^, mag), as described 
in the text. 



initions are taken from Luppino & Kaiser (1997). Note 
that we approximate the matrix P 7 by a scalar equal to 
half its trace. Since measurements of Tr P 7 from indi- 
vidual galaxies are noisy, we follow Fu et al. (2008) and 
fit it as a function of galaxy size and magnitude, which 
are more robustly observable galaxy properties. 

Following Hoekstra et al. (2000), we weight the shear 
contribution from each galaxy as 



p2 



(2) 



where a e ,i is the error in an individual cllipticity measure- 
ment obtained via the formula in Appendix A of Hoek- 
stra et al. (2000), and Co ~ 0.278 is the dispersion in 
galaxies' intrinsic ellipticities. 

3.4. Calibration of Multiplicative Shear Measurement 

Biases 

We exploit the opportunity that the CFHTLS-Deep D2 
field includes the HST COSMOS survey field, and ver- 
ify the calibration of our shear measurement pipeline for 
ground-based data against an independent analysis of the 
much higher resolution space-based data (Leauthaud et 
al. 2007, 2010). We stack subsets of the CFHTLS-Deep 
D2 imaging to the same depth as the CFHTLS-Wide 
survey and analyze it using the same pipeline applied 
to the CFHTLS-Wide Wl field. Since any galaxy seen 



IQ=0.90, m — 0.21±4.61E-3. c — 1 .33E-2±3.39E-3 
IQ-0.69. m--0.21±3.47E-3. c=- 1.09E-2±2.58E-3 
IQ=0.57. m = -0.21±3.02E-3. c=-1.05E-2±2.25E-3 
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Figure 8. Robustness of the calibration of our shear measurement 
as a function of image quality. Linearly-spaced contours compare 
our shear measurements of galaxies in subsets of CFHTLS-Deep 
imaging with varying mean seeing (black solid: 0".90, red dashed: 
0".69, blue dotted: 0".57) to measurements from the Hubble Space 
Telescope. Shears are consistently underestimated by our pipeline, 
but the calibration is remarkably robust. 



by CFHT is very well resolved by HST, and imaged to 
very high signal-to-noise ratio by the COSMOS survey, 
the space-based shear measurements require only negli- 
gible PSF correction and suffer from only negligible shot 
noise. A consistent shear measurement between ground 
and space for this subset of galaxies would therefore in- 
dicate a robust shear measurement across the CFHTLS. 

Multiplicative shear measurement biases m are the 
most problematic for circularly-symmetric cluster mea- 
surements. Multiplicative biases cannot be internally 
diagnosed within a shear catalog, so our comparison 
against external data is most useful for checking that 
m is sufficiently small that it corresponds to a bias 
smaller than our statistical errors. Within the KSB 
framework, difficulties in shear calibration mainly rest 
in measurement of the shear polarizability P 7 , so we 
first investigate different possibilities for fitting P 7 across 
a galaxy population. Figure [7] compares shear mea- 
surements from a subset of the CFHTLS-Deep imag- 
ing with mean seeing 0".69 (similar to the mean see- 
ing in our survey) stacked to the depth of the CFHTLS- 
Wide imaging against shear measurements obtained from 
HST. The dashed lines show the best-fit linear relations 
^,cfht _ ^ _|_ m ^HST _j_ c ^ which are obtained using 

a total least-squares fitting method (e.g. Kasliwal et al. 
2008) that accounts for the noise present in both shear 
catalogs. The top- left panel shows the CFHT shear 
measurements with P 7 naively obtained from each raw, 
noisy galaxy without any fitting. This results in a large 
bias on shear measurements and a large amount of ex- 
tra noise. The top-right panel shows the shear mea- 
surements if P 7 is fitted as a function of galaxy size, 
P.y(r h ) = a Q + air h - 



a,ir\. The bottom-left panel shows 
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Figure 9. Mean shear measurements from galaxies in i'-band ob- 
servations of the entire CFHTLS-Wide Wl field. In the absence 
of additive systematics, these should be consistent with zero. In 
practice, they always remain within the dashed lines than indicate 
an order of magnitude lower than the 1%-10% shear signal around 
clusters. Upper and lower panels show components 71 and 72, re- 
spectively. Left, middle, and right panels show trends as a function 
of galaxy magnitude, size, and detection signal-to-noise. 

shears if P 7 (r^, mag) = ao+axrh,+a2rf l +a^m v (Fu et al. 
2008). The bottom-right panel shows the matched shear 
with P 7 (th , mag) the best-fit rational function 



a + a\m v + a 2 mf, + a 3 r h 
1 + a 4 m, y + a 5 ml + a&rh + a^r\ 



(3) 



In this example, the coefficients are ao = 25.07, a± = 
-2.19, a 2 = 0.045, a 3 = 0.53, a 4 = 0.58, a 5 = -0.022, 
a 6 = —0.85, and 07 = 0.14. The more sophisticated fits 
produce a shear catalog that is a marginally better match 
to the reliable HST measurements, and this is even more 
true if we redo the analysis using the full CFHTLS-Deep 
depth, in which galaxies are fainter and smaller. We 
henceforth choose to adopt the rational function fit to 
Pj for all subsequent analyses, obtaining new best-fit 
coefficients for each pointing. 

To quantify the performance of our shear measurement 
pipeline as a function of image quality, we stack subsets 
of the CFHTLS-Deep D2 imaging with low-, medium-, 
and high-seeing to the same depth as the CFHTLS-Wide 
survey, and analyze each separately (Figured]). We find 
that our CFHTLS pipeline consistently underestimates 
shear, but that the calibration is remarkably robust to 
seeing conditions. We can therefore simply recalibrate 
our pipeline for all images by multiplying all measured 
shears by 1/(1 - 0.21). 

3.5. Assessment of Residual Additive Shear Systematics 

Additive shear measurement systematics c generally 
cancel out in circularly-symmetric cluster measurements 
(Mandelbaum et al. 2006). However, to double-check 
for significant additive systematics, we first measure the 
mean shears (7) across all 72 pointings of the CFHTLS- 
Wide Wl field. Figure [9] demonstrates that the mean 
shear is consistent with zero as expected, for galaxies of 
all sizes, magnitudes, and signal-to-noise ratios. 

We also look for residual systematics left in the weak 
lensing cosmic shear signal due to imperfect PSF correc- 
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Figure 10. The cross-correlation between shear measurements 
and stellar ellipticities, as a function of the separation between 
galaxies and stars, averaged throughout the CFHTLS-Wide Wl 
field. If all residual influence of the observational PSF has been 
successfully removed from the galaxy shape measurements, the red 
(lower) points should be consistent with zero. 



tion. Figure [10] shows the correlation £ sys between the 
corrected shapes of galaxies and the uncorrected shapes 
of stars. Following Bacon et al. (2003) and Massey et 
al. (2005), we normalize the star-galaxy ellipticity corre- 
lation by the uncorrected star-star ellipticity correlation 
to assess its impact on shear measurements 



£sys(#) 



< e*(x) 7 (x + 6>) > 2 

< e*(x) e*(x + 0) > 



(4) 



where e* is the ellipticity of the stars before PSF cor- 
rection and 7 is the shear estimate from galaxies. We 
find that our PSF correction is well within requirements 
for our analysis because on cluster scales 1-5 arcmin, 
the amplitude of £ S ys is at least one order of magnitude 
smaller than the cosmic shear signal 



£± = &t(0)±&*(0) 



1 

2^ 



eP K (£)J M9)6£, (5) 



where £tt(#) (£xx(#)) are the correlation functions be- 
tween components of shear rotated tangentially (at 45°) 
to the line between pairs of galaxies separated by an an- 
gle 8 and Jo, J4 are Bessel functions of the first kind. 

Because gravitational lensing is achromatic while sys- 
tematics are typically not, we can also assess the robust- 
ness of our measurements by comparing shears measured 
from independent imaging acquired in multiple bands. 
The first attempt at comparing multicolor shear mea- 
surements was made by Kaiser et al. (2000) using the 
CFHT12K camera. The I and V bands showed signif- 
icantly different signals that were inconsistent with the 
change in redshift distribution between the two filters. 
After a great deal of algorithmic progress, Semboloni et 
al. (2006) obtained consistent shear measurements from 
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Figure 11. Gravitational lensing is achromatic, so measurements 
of galaxy shapes from imaging in different colors should on average 
be consistent. This shows a comparison of shear measurements 
obtained from CFHTLS-Wide r'-band and i'-band imaging of the 
whole Wl field. The dashed line shows the best-fit linear relation 
7* = (in — 1)7[ + c. 



i'-band and r'-band CFHTLS-Deep data. Gavazzi & 
Soucail (2007) extracted consistent shear from the g', r', 
z' and i' bands of CFHTLS-Deep. Gavazzi et al. (2009) 
also measured consistent values of the PSF-corrected el- 
lipticities of central Coma cluster galaxies in MegaCam- 
u* and CFH12k-7 bands. Our analysis pipeline mea- 
sures shears in the independent CFHTLS r'-band and 
i'-band imaging that are consistent within the ~ 0.1 
rms noise (Figure ITT]) . To maximize the total number 
of galaxies in our shear catalog, we therefore combine 
r'-band and i'-band measurements. Only the unique r'- 
band galaxies are added to the i'-band catalog. The 
combined catalog includes ~ 2.66 million galaxies, with 
n g ~ 14.5 arcmin -2 for lensing measurements. 

We shall continue testing for systematics at each stage 
of our analysis, by checking that the shear signal behaves 
as expected, and is consistent with external data sets. An 
important example of this is the (nonphysical) B-mode 
signal, which we shall compute wherever we reconstruct 
the (physical) E-mode. Pure gravitational fields produce 
zero £?-mode for isolated clusters and only tiny B-modes 
through coupling between multiple systems along adja- 
cent lines of sight (Schneider et al. 2002). The -B-mode 
signal corresponds to the imaginary component of P K (l)\ 
it can be conveniently measured by rotating all galaxy 
shears through 45° then remeasuring the -E-mode signal 
(Crittenden et al. 2002). 

4. MASS RECONSTRUCTIONS 

4.1. Kaiser-Squires Inversion and Masking 

The shear field Ji(6) is sparsely and noisily sampled 
by measurements of the shapes of galaxies at positions 9. 
The smooth, underlying shear field 7i(0) can be written 
in terms of the lensing potential 4>{0) as 



li = \{dl-dt)<i>, 
72 = <9i<9 2 0, 



(6) 
(7) 



where the partial derivatives di are with respect to 
The convergence field k(0), which is proportional to the 
mass projected along a line of sight, can also be expressed 
in terms of the lensing potential as 



n=\{dl+dl) 



(8) 



We shall reconstruct the convergence field from our 
shear measurements via the Kaiser & Squires (1993) 
(KS93) method. This is obtained by inverting Equa- 
tions ([6]) and (JJJ in Fourier space: 



7» 



= Pik 



(9) 



for i = 1,2, where the hat symbol denotes Fourier trans- 
forms, we define k 2 



k\ 



fc| and 



A (*0 



A(fc) 



kl 



k 2 



k 2 
2kik 2 



k 2 



(10) 



(11) 



This inversion is non-local, so we deal with masked re- 
gions of the shear field by masking out the same area in 
the convergence field, plus a 1'.5 border. We shall ignore 
any signal within these regions, and set the convergence 
to zero in relevant figures. 

For the finite density of source galaxies resolved by 
CFHT, the scatter of their intrinsic ellipticities means 
that a raw, unsmoothed convergence map k(9) will be 
very noisy. Following Miyazaki et al. (2002), we smooth 
the convergence map by convolving it (while still in 
Fourier space) with a Gaussian window function, 



(12) 



As shown by van Waerbeke (2000), if different galax- 
ies' intrinsic ellipticities are uncorrelated, the statistical 
properties of the resulting noise field can be described by 
Gaussian random field theory (Bardeen et al. 1986; Bond 
& Efstathiou 1987) on scales where the discreteness ef- 
fect of source galaxies can be ignored. The Gaussian field 
is uniquely specified by the variance of the noise, which 
is in turn controlled by the number of galaxies within a 
smoothing aperture (Kaiser & Squires 1993; Van Waer- 
beke 2000) 

o-l 1 

(13) 



-noise 2 27r 02 rls > 

where a e is the rms amplitude of the intrinsic elliptic- 
ity distribution and n g is the density of source galaxies. 
We define the signal-to-noise ratio for weak lensing de- 
tections by 

v=-^-. (14) 

f^noisc 

To define the noise level in theoretical calculations of v, 
we adopt a constant effective density of galaxies equal 
to the mean within our survey. For observational cal- 
culations of v, we use the mean galaxy density in each 
pointing — but do not consider the non-uniformity of the 
density within each field due to masks or galaxy cluster- 
ing. 
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Figure 12. Distribution of foreground mass in the Wl+2+3 pointing, reconstructed from shear measurements via the KS93 method. 
Left: the physical B-mode convergence signal. Right: the B-mode systematics signal, created by rotating the shears by 45° then remaking 
the map. Contours are drawn at detection significances of 3<r, 4er, and 5<r, with dashed lines for negative values. 



It turns out that a simple Gaussian filter of width 
Oq k, V is close to the optimal linear filter for cluster 
detection, and this choice has been extensively studied 
in simulations (White et al. 2002; Hamana et al. 2004; 
Tang & Fan 2005). Because of our relatively low source 
galaxy density, the galaxies' random intrinsic shapes will 
produce spurious noise peaks, degrading the complete- 
ness and purity of our cluster detection. To reduce con- 
tamination, we repeat our mass reconstruction using two 
smoothing scales 6 a = 1' and 0q = 2'. The map with 
greater smoothing will be less noisy; to help remove spu- 
rious peaks from the higher resolution map, we consider 
only those peaks detected above a signal-to-noise thresh- 
old in both maps. 

Figure [T2] shows the reconstructed convergence field 
corresponding to foreground mass in the Wl+2+3 point- 
ing. The left panel shows the B-mode reconstruction 
with KS93 method after smoothing by a larcmin Gaus- 
sian kernel. This contains several high signal-to-noise ra- 
tio peaks, while the associated B-mode systematics mea- 
surement in the right panel is statistically consistent with 
zero, with fewer peaks. As weak lensing produces only 
curl-free or B-mode distortions, a detection (significant 
above statistical noise) of curl or B-mode signal would 
have indicated contamination from residual systematics, 
e.g. imperfect PSF correction. 

4.2. Large-scale Lensing Mass Map 
A reconstructed "dark matter mass" convergence map 
for the entire 64 deg 2 CFHTLS-Wide Wl field is pre- 
sented in Figure [T51 We detect 301 peaks with v > 3.5 
in maps with both smoothing scales 9q = 1' and 0q = 2'. 
The same information is reproduced in Figure [T4l with 
9q = 6' and after multi-scale entropy restoration filtering 



Table 2 

The Number of Local Maxima and Minima in the Convergence 
Map of the CFHTLS-Wide Wl Field, as a Function of Smoothing 

Scale. 



Smoothing 
Scale 8q 


B-mode 
v > 3.5 


B-mode 
v < -3.5 


B-mode 
v > 3.5 


B-modc 
v < -3.5 


0'.5 


1512 


1270 


1244 


1033 


l'.O 


543 


445 


361 


282 


2'.0 


281 


233 


148 


126 



(MRLens; Starck et al. 2006), to better display large- 
scale features. The MRLens filtering effectively sup- 
presses noise peaks, but results in non-Gaussian noise 
that complicates the peak selection (Jiao et al. 2011), so 
we shall not use it further. 

To assess the reliability of this map, we shall first inves- 
tigate the statistical properties of local maxima and min- 
ima. Figure [15] shows the distribution of peak heights, 
as a function of detection signal-to-noise. The bimodal 
distribution in both the B-mode and B-mode signals is 
dominated by positive and negative noise fluctuations, 
but an asymmetric excess in the B-mode signal is appar- 
ent at both v > 3.5 and, at lower significance, v < —3.5. 
The amplitude, slope and non-Gaussianity of this excess 
are all powerful discriminators between values of param- 
eters in cosmological models (Pires et al. 2009). Positive 
peaks correspond mainly to dark matter halos around 
galaxy clusters. Local minima could correspond to voids 
(Jain & van Waerbeke 2000; Miyazaki et al. 2002), but 
the large size of voids is ill-matched to our 9 a — V filter 
width, and their density contrast can never be greater 
than unity, so this aspect of our data is likely just noise. 

Figure [16] recasts the peak distribution into a cumu- 
lative density of positive maxima or negative minima. 
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Figure 13. Reconstructed "dark matter mass" convergence map for the entire 64 deg 2 CFHTLS-Wide Wl field. This has been smoothed 
by a Gaussian filter of width 6q = 1'. Black contours are drawn at detection signal-to-noise ratios v = 3.0, 3.5, 4.0, and red contours 
continue this sequence from v > 4.5. 



As expected, we find a non-Gaussian mass distribution 
with more highly significant positive maxima (corre- 
sponding to mass overdensities) than highly significant 
negative minima (see Table [2]) . Analytic predictions of 
peak counts are also overlaid. Following van Waerbeke 
(2000), dot-dashed lines show the expected density of 
pure noise peaks, and dotted lines show the expected 
number of true dark matter halos. Predictions from Fan 
et al. (2010), which also take into account the effect of 
noise on the heights of true peaks and the clustering of 
noise peaks near dark matter halos, are shown as cir- 



cles with error bars. In these theoretical calculations, 
we model the population of background galaxies as hav- 
ing an intrinsic ellipticity dispersion a e — 0.278, density 
n g = 14.5 arcmin -2 and the redshift distribution from Fu 
et al. (2008) in a ACDM universe. At v > 4.5, it appears 
that theory may begin to predict more peaks than are 
observed. However, these are very small number statis- 
tics, and our observations are consistent with analytical 
predictions within Poisson noise. At very low v, the num- 
ber of peaks is washed out and actually decreases when 
noise is superimposed, because it is impossible to extract 
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Figur e 14 . Reconstructed "dark matter mass" convergence map for the entire 64 deg 2 CFHTLS-Wide Wl field showing the same data as 
Figure IT3l but smoothed with a Gaussian filter of width 8q = 6' (left) and multi-scale wavelet filtering (right) to highlight the large-scale 
features. 



very low-v peaks. However, if larger and deeper surveys 
still find fewer low-is or high-^ peaks than expected, and 
systematic effects such as the consequences of masked 
regions are more fully understood, it may indicate a cos- 
mology with, e.g., lower VL m and erg than the values used 
for our predictions. 

Our main conclusions about the distribution of conver- 
gence peaks are as follows. 

• Peak counts detected in CFHTLS-Wide are consis- 
tent with predictions from a ACDM cosmological 
model, once noise effects are properly included (van 
Waerbeke 2000; Fan et al. 2010). 

• The convergence field is non-Gaussian. The excess 
of local maxima with v > 3.5 compared to local 
minima with v < —3.5 is also consistent with mod- 
els (Miyazaki et al. 2002; Gavazzi & Soucail 2007). 

• Noise peaks dominate the expected peak counts 
due to cosmological weak lensing below ;$ 3. 

We expect that weak lensing peak counts will become 
reliably employed to constrain cosmological parameters 
in future lensing analyses. Although it is difficult to 
remove the contributions of noise from intrinsic galaxy 
shapes and the projection of large scale structures, these 
effects can be analytically predicted. Pushing these pre- 
dictions into the low signal-to-noise regime might also 
allow constraints to use a much higher number density 
of less massive peaks, tightening predictions by reducing 
Poisson noise. 

5. OPTICAL/X-RAY COUNTERPARTS 

We shall now compare our weak lensing peak detec- 
tions with catalogs of overdensities detected via optical or 



X-ray emission. We adopt the K2 optical cluster catalog 
constructed using photometric redshifts from the same 
CFHTLS-Wide Wl imaging (Thanjavur et al. 2009), and 
the XMM-LSS X-ray selected cluster catalog (Adami et 
al. 2011), which partially overlaps our survey field. 

We define the purity f p of our blind weak lensing clus- 
ter search as the fraction of peaks above a given detec- 
tion threshold i>th that are associated with an optically 
detected cluster 

f _ A^matchcdC^ > ^th) -i 

h ~ N{v>v th ) ■ [lb) 

The expected purity depends upon the survey depth 
(galaxy density), systematics, and the extent of multi- 
band or spectroscopic follow-up. Using numerical simu- 
lations with n g ai = 30 arcmin -2 , Hamana et al. (2004) 
predict a purity of more than 60% for convergence peaks 
with v > 4. The Bonn Lensing, Optical, and X-ray se- 
lected galaxy clusters (BLOX) simulations by Dietrich et 
al. (2007) with n ga i = 18 arcmin -2 also predict that 75% 
of matches between convergence map peaks and massive 
halos are within 2'. 15. In practice, Miyazaki et al. (2007) 
achieve 80% purity for the ~ 100 peaks in Subaru con- 
vergence maps with v > 3.69. Gavazzi & Soucail (2007) 
obtain ~ 65% purity from 14 peaks in CFHTLS-Deep 
with v > 3.5. Schirmer et al. (2007) obtain ~ 45% purity 
for the 158 possible mass concentrations identified in the 
Garching-Bonn Deep Survey (GaBoDS) at v > 4, con- 
sistent with an earlier evaluation of a subsample of the 
survey (Maturi et al. 2007). Geller et al. (2010) find only 
~ 33% purity by combining the Deep Lens Survey with 
the Smithsonian Hectospec Lensing Survey (SHELS) (see 
Table [2]) . 

Figure [17] shows a subset of our cluster search in the 
representative Wl+2+3 pointing, which also overlaps 
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Table 3 

The Detection Purity of Various Weak Lensing Cluster Surveys Described in the Literature. 



Weak Lensing Cluster Sample Comparison Sample Purity 



Hamana et al. (2004) 


n gal = 30 arcmin 2 , v > 4 


Halos (simulations) 


60% 


Dietrich et al. (2007) 


n gal = 18 arcmin -2 , 


Halos (simulations) 


75% 


Miyazaki et al. (2007) 


v > 3.69 


X-ray clusters (XMM-LSS) 


80% 


Gavazzi & Soucail (2007) 


"gal ~ 35 arcmin -2 , v > 3.5 


Photometric clusters 


~ 65% 


Schirmer et al. (2007) 


v > 4 


BCG 


~45% 


Geller et al. (2010) 


v > 3.5 


Spectroscopic clusters (SHELS) 


- 33% 



10' 




10' 



1000 7 




Figure 15. Numbers of local maxima (solid line) and minima 
(dashed line) in our E-mode (black) and B-mode (red) convergence 
map of the CFHTLS-Wide Wl field, with smoothing scales B G = V 
(top-panel) and 8q = 2' (bottom-panel). Local maxima can still 
have a slightly negative peak height if they occur along the same 
line of sight as a negative noise fluctuation (or a huge void), and 
local minima can similarly have a slightly positive peak height. To 
eliminate spurious noise peaks, we shall mainly consider maxima 
or minima with \u\ > 3.5. In this regime, there is an excess of local 
maxima over local minima, and an excess of E-mode peaks over 
B-mode peaks (see Table [2J. 
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Figure 16. Cumulative density of local maxima 7V(> v) (solid 
line) and corresponding density of local minima N(< —v) (dashed 
line) in the CFHTLS-Wide Wl B-mode (black) and B-mode (red) 
convergence map with smoothing scale 9q = 1' (top-panel) and 
8q = 2' (bottom-panel). Error bars are simply la uncertainties 
assuming Poisson shot noise. The dot-dashed line shows the pre- 
diction from Gaussian random field theory (van Waerbeke 2000). 
Circles with additional error bars and the dotted curves show an 
analytical prediction in a ACDM universe (Fan et al. 2010), with 
and without the influence of random noise. The excess of positive 
maxima over negative minima demonstrates the non-Gaussianity 
of the convergence field. 



with the XMM-LSS survey (Pacaud et al. 2007; Adami et 
al. 2011). The observed weak lensing peak positions may 
not coincide exactly with cluster centers defined from 
optical emission, because of a combination of noise, sub- 
structure, and physical processes associated with cluster 
mergers (Fan et al. 2010; Hamana et al. 2004). We there- 
fore search for matched pair candidates in K2 within a 



3'.0 radius of peaks that appear in both the Bq = 1' and 
6q = 2' lensing map. This search radius is chosen to 
be larger than the smoothing scale, but smaller than the 
angular virial radius of a massive cluster at0.1<z<0.9 
(Hamana et al. 2004). If more than one pair exists within 
3'.0, we adopt the closest match as the primary candi- 
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Figure 17. Reconstructed lensing convergence signal-to noise 
map for the Wl+2+3 pointing, plus overlays showing optically 
and X-ray selected cluster counterparts. The smoothing scale of 
the background map is Qq = 1' '. Symbols indicate the positions of 
+: lensing peaks detected with v > 3.5 in the Qq = 1' map, 
□: lensing peaks detected with v > 3.5 in the 9q = 2' map, 
A: optically-detected clusters in the K2 catalog, and 
X: X-ray selected clusters found in XMM-LSS (Adami et al. 2011). 



Table 4 

Cluster Search Purity as a Function of Peak Height Threshold 



f-th 


N{u > i/ th ) 


^matched (f > ^th,K2) 


f P 


3.5 


301 


126 


42% 


4.0 


125 


67 


54% 


4.5 


51 


30 


59% 



date. 

We obtain 126 matches between weak lensing peaks 
with v > 3.5 and optical K2 clusters: corresponding to 
a purity of 42%. This is lower than in CFHTLS-Deep 
(Gavazzi & Soucail 2007) because of the much lower 
source galaxy density. The purity is listed as a function of 
detection threshold t'th in Table HI and the complete cat- 
alog of matched pairs is presented in Table [5j Figure [18] 
shows the separations between the matched weak lensing 
peaks and K2 centers. The left panel shows the offsets 
for peaks with v > in maps with various smoothing 
scales; the separation between peaks typically matches 
the smoothing scale. The right panel shows the offsets 
for only the reliably-detected peaks with high v : the fi- 
nite number of clusters is too low to draw solid conclu- 
sions about the typical separation. Figure [T^l shows the 
redshift distribution of matched clusters, relative to the 
overall K2 sample. Our lensing selection preferentially 



detects clusters at 0.2 < z < 0.4, and becomes inefficient 
above z > 0.5. 

We also compare our v > 3.5 lensing peaks with X- 
ray observations from the XMM-LSS survey. Adami et 
al. (2011) present 66 spectroscopically confirmed clusters 
(0.05 < z < 1.5) within the 6 deg 2 XMM-LSS survey. 
This partially overlaps with the CFHTLS-W1 field (53 
X-ray clusters are within Wl). In this overlap region, we 
find 31 lensing peaks, 7 of which are within 5'.0 of X-ray 
clusters. Note that we increase the distance threshold for 
matches because of additional noise in the X-ray centers 
and the common phenomenon of separation between the 
gravitational field and X-ray gas — Shan et al. (2010) 
found that 45% of a sample of 38 clusters had X-ray 
offsets > 10". Indeed, in our new sample, we also find 
that the offsets between weak lensing and X-ray centers 
are always comparable to or much larger than the offset 
between weak lensing and optical centers. 

Miyazaki et al. (2007) also perform a weak lensing 
analysis of 0.5 deg 2 of the XMM-LSS survey. They find 
15 lensing peaks with optical counterparts, of which 10 
match X-ray selected clusters (Adami et al. 2011). Many 
of the X-ray clusters are simply at too high a redshift 
to be detected by the CFHTLS lensing data. However, 
three of these clusters are detected in both our analysis 
(c77, c92 and c93 in Table [6]) and that of Miyazaki et al. 
(2007). 

6. TOMOGRAPHIC ANALYSIS OF LENSING PEAKS 

Since multi-band photometric redshift estimates are 
available for 76% of the source galaxies, we shall now per- 
form a 3D, tomographic shear analysis around all v > 3.5 
cluster candidates with an optical or X-ray counterpart. 
This process yields an estimate of the cluster redshift 
(and mass) independently of its visible emission. It also 
further cleans the cluster catalog of spurious peaks cre- 
ated by either noise or the projection of multiple small 
systems along a single line of sight. 

We fit the 3D shear signal around each cluster can- 
didate to both a singular isothermal sphere (SIS) and 
a Navarro et al. (1996; NFW) model, with the cluster 
lens as an additional free parameter. We assume that 
the center of each cluster candidate is at the position 
of the brightest central galaxy (BCG). This may not be 
optimal (Johnston et al. 2007), but it is much more pre- 
cisely known than the peak of the lensing signal (Fan et 
al. 2010). For each source galaxy, we adopt the best- 
fit redshift from the cleaned Arnouts et al. (2010) cata- 
log. Following Gavazzi & Soucail (2007), we measure the 
mean tangential shear in a l'-5' annulus from the clus- 
ter center, excluding the core to minimize dilution of the 
signal from any cluster member galaxies with incorrect 
photometric redshifts. 

For an SIS model, the component of shear tangential 
to the cluster center is 



7t(0, z s ;z t ) 



Dis 4w 2 1' 



29 



(16) 



where 9 is the angular distance from the center. We 
simultaneously fit the unknown lens redshift and charac- 
teristic cluster velocity dispersion a v by minimizing 



xli S (zi,a v ) = 



E 



(it.i - 7SIS,i(zz, cr tomo )) 



(17) 
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Figure 18. Distribution of the offsets between matched pairs of weak lensing peaks and K2-detected clusters. Left: offsets for all peaks 
in maps with smoothing scales of 0'.5 (solid black), 1' (dotted red), and 2' (dashed blue). Right: offsets for peaks with v > 3 (solid black 
histograms), v > 3.5 (dotted red histograms) and v > 4 (dashed blue histograms), all with the smoothing scale 8 = 1 arcmin. 
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Figure 19. Cluster redshift distribution for the matched clusters 
of weak lensing and K2 (black solid histogram) and the total K2 
detected clusters with K2 detection significance (r — i) > 3 (red 
dashed histogram). 

where Wi = w(zi, z s ^) and <j 1 ei is given by Equation. (2). 

The full NFW model has two free parameters, but we 
assume the Bullock et al. (2001) relation between con- 
centration c and virial mass M v j r seen in numerical sim- 
ulations 

c / M ■ \-° 13 
c NFW (M vir , z) = — [y^J^ ) , (18) 

where c* = 8. As for the SIS model, there is then only 
one free parameter, and we fit the shear field for lens 
redshift and halo mass by minimizing 

2 / , f \ (7M -7NFw(^,M vir )) 2 

X N Fw( Z *> M vir) = 2^ 2 — • ( 19 ) 

i 

To aid comparison with other work, for each cluster we 
calculate both the virial mass M v ; r and the mass M200 



enclosed within a radius r2oo in which the mean density 
of the halo is 200 times the critical density at the redshift 
of the cluster. 

Figurel20l illustrates our 3D tomographic results on one 
cluster (identification c3). The mean tangential shear 
is consistent with zero for z s < 0.26. The subsequent 
increase with redshift is clear and allows for an unam- 
biguous identification of the lens redshift. Error bars are 
derived from the scatter in observed ellipticities (intrin- 
sic+measurement error) , as determined by Equation (TIJ . 
The observed tangential shear profile is consistent with 
either an SIS or NFW model. Due to our low density of 
source galaxies and our exclusion of galaxies within the 
central 1', we have insufficient signal-to-noise for individ- 
ual clusters to distinguish between the two models, which 
differ most noticeably near the core. The amplitude of 
the systematic £?-mode signal, computed by rotating all 
shear estimates by 45°, is always at least one order of 
magnitude smaller than the tangential shear, and it os- 
cillates about zero. 

Figure [21] shows the signal-to-noise contours of the con- 
vergence signal reconstructed around four clusters from 
the 3D shear signal. 

Our tomographic analysis confirms the identification of 
85 clusters with x^cd = X 2 /dof < 3. The mean redshift 
and velocity dispersion of these clusters are (z c ) — 0.36 
and (a v ) — 658.8km s _1 . Their full properties are listed 
in Table [SJ Reassuringly, we find that the inferred lens 
redshifts are effectively identical for either SIS or NFW 
profile fits, and are consistent (although noisy) with the 
(independent) photometric redshifts of the cluster mem- 
ber galaxies. Figure 1221 compares the redshift estimates 
from tomographic gravitational lensing and optical spec- 
troscopy of member galaxies for these 85 clusters. 

At the faint end of our source galaxy sample, photo- 
z estimates will be unreliable because of noise in the 
photometry and degeneracies in the broad-band colors 
of galaxies with different spectral energy distributions 
at different redshifts. This will show up as a "double 
peak" in the posterior probability of the redshift dis- 
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Figure 20. Results of a tomographic weak lensing analysis around 
one example peak (c3 in our catalog), which has an optical coun- 
terpart at redshift z = 0.28. The solid and dashed lines are the 
best-fit SIS and NFW models. Top: projection in the redshift di- 
rection showing the characteristic increase with redshift of a real 
signal. The best-fit lens redshift is 0.26^q 'J for an SIS model 
and 0.26^q ijj f° r an NFW model. Middle: radial profile projected 
onto the plane of the sky. Black circles show the B-mode tangential 
shear signal. Red circles (triangles) show positive (negative) values 
of the B-mode systematic signal, which oscillates about zero. Bot- 
tom: joint redshift-mass constraints from the best fit NFW model. 
Contours show 68%, 95% and 99% confidence limits. 



tribution. Using only the best-fit peak might be ran- 
domly picking whichever of these peaks is higher be- 
cause of noise. This often biases lensing analyses be- 
cause the expected lensing signal may be much higher at 
one redshift than the other. To check for such an effect, 
we redo the tomographic analysis without any photo- 
zs that have a double peak (Arnouts et al. 2010). For 
cluster c3 in our catalog, we get very similar fit results: 
zsis = 0.26tg;^ and a v = 724.0t^J for the SIS model, 
and z NFW = 0.26to;™ and m NFW = 4.141°;^ for the 
NFW model. This suggests that the fit is not very sensi- 
tive to that population of galaxies with "double peaks" 
in photometric redshifts. 

In cases where the data are poorly fit (x 2 od > 3) by 
a 3D lensing signal, the inferred velocity dispersion or 
virial mass are typically very small (often consistent with 
zero) . These systems are probably spurious peaks due to 
noise or projection effects. In addition, two peaks (c71 
and cl04) have an unphysically high velocity dispersion 
ctomo > 1300 km s _1 . After careful examination of the 
image data, we found that they lie near two strongly 
saturated stars whose flux extends beyond the masked 
regions, possibly degrading galaxy shape measurements. 

One important goal of cluster lensing is to measure the 
total mass of systems, which can also be estimated from 
the velocity dispersion of its member galaxies or from 
its X-ray emission, under various assumptions about the 
state of the intra-cluster medium and hydrostatic equi- 
librium (e.g., Bahcall et al. 1995; Carlberg et al. 1996; 
Carlstrom et al. 2002). 

Figure [23] shows the relation between one mass proxy 
obtained from X-ray observations and another obtained 
from weak lensing, for the 4 of 7 clusters in our sample 
with an X-ray counterpart and an SIS velocity disper- 
sion parameter well constrained by lensing. These are 
overlaid on the results of another lensing-selected cluster 
sample (Hamana et al. 2009) and two X-ray selected clus- 
ter samples (Cypriano et al. 2004; Hoekstra 2007). Even 
combining these samples, the scatter is large. However, 
the consistency of the scatter between the four samples 
suggests that neither detection method produces a strong 
selection bias. Note that, since the density (and shear) 
profile of real cluster is not necessarily a single power 
law, the best-fit SIS model may depend on details of the 
fitting method, and the range over which data are fit. 
Our results therefore are somewhat method-dependent. 
A corollary of this issue is that it might also be possible 
to minimize scatter in the a xray -at om o relation by opti- 
mizing the tomographic lens fitting method. 

7. CONCLUSIONS 

We have presented a tomographic weak gravitational 
lensing analysis of the completed 64 deg 2 CFHTLS-Wide 
Wl field, demonstrating some of the power of lensing to 
probe mass in galaxy clusters. We measured the shapes 
of distant galaxies using the KSB shape measurement 
method, which we verified against shape measurements 
from high-resolution HST imaging of an overlapping sky 
area. We also obtained consistent shape measurements 
using two independent imaging bands. The level of resid- 
ual shape measurement systematics is an order of mag- 
nitude lower than the 1%-10% shear signal expected in 
galaxy clusters, so this analysis is acceptable for cluster 
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Figure 21. Composite 3x3 arcmin 2 CFHTLS g' , r' , i' color images for four clusters detected in both weak lensing and the K2 optical 
catalog. Contours show signal-to-noise in reconstructed convergence, starting at v = 3.0 and increasing in steps of 0.5 and the images are 
centered on the convergence peak. Top-left: cluster c3 with v = 5.395. Top-right: cluster c49 with v = 4.728. Bottom-left: cluster c21 
with v = 4.142. Bottom-right: cluster cl03 with v = 3.725. Candidate c21 is included here as an example of a 2D lensing peak that is 
probably spurious: it does not appear obviously associated with an overdensity of galaxies in the optical imaging, and is not well-fit by a 
3D shear signal. 



studies. 

We have reconstructed the largest contiguous "dark 
matter map" convergence field to date, using two dif- 
ferent smoothing scales to help remove spurious noise 
peaks. From this map, we performed the largest lensing- 
selected blind cluster search to date, finding 301 local 
maxima in the lensing map with detection significance 
v > 3.5. Once sources of noise are properly modeled 
from the intrinsic shapes of galaxies and substructure, 
this is consistent with predictions from a ACDM cos- 
mology. Note that our theoretical calculations (Fan et 
al. 2010) do not consider projection effects of structures 
along the line of sight. 



To identify counterparts of the weak lensing peaks, 
we match our cluster candidates to the K2 galaxy clus- 
ter catalog, created using photometric redshift estimates 
across the CFHTLS- Wide Wl field (Thanjavur et al. 
2009). Of the 301 peaks with v > 3.5, 126 have a cor- 
responding optically-detected BCG within 3'.0. In the 
(much smaller) survey area that overlaps the XMM-LSS 
survey, we also find matches for seven lensing peaks with 
X-ray selected clusters. Thus, many of the candidate 
peaks are indeed likely just noise. 

Tomographic weak lensing techniques dramatically im- 
prove standard 2D algorithms. In a full 3D lensing anal- 
ysis of the v > 3.5 peaks, we further distinguish real clus- 
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Figure 22. Comparison between cluster redshifts derived directly 
from the photometric redshifts of cluster member galaxies z op t and 
from weak lensing tomography ztomo- This figure includes all 85 
v > 3.5 weak lensing peaks with optical counterparts, whose 3D 
shear signal is consistent at Xred ^ ^ with the expected increase as 
a function of redshift. 
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Figure 23. Scaling relation between X-ray and weak lensing 
tomographic measurements of the velocity dispersion of clusters. 
Black filled circles show our data for four lensing-selected clusters. 
Red filled squares show the lensing-selected clusters of Hamana et 
al. (2009). Blue open circles and green open squares denote the 
X-ray selected clusters from Cypriano et al. (2004) and Hockstra 
(2007) respectively. 



ters from noise fluctuations, and confirm (at Xred < 3) 
the identification of 85 clusters. Importantly, we ob- 
tain independent measurements of the cluster lens red- 
shifts, which are consistent with the redshifts of their 
previously-identified optical counterparts. For each clus- 
ter, we fit NFW and SIS radial profiles to the lensing 
data to measure the mass or velocity dispersion <7t om o- 
The clusters' mean redshift and velocity dispersion is 
(z c ) = 0.36 and (er c ) = 658.8 km s _1 . Weak lensing 
measurements of the total mass in the four of our clusters 
with X-ray counterparts are also in reasonable agreement 
with mass estimates obtained from X-ray emission. Fu- 
ture surveys, such as DES, LSST, KDUST and EUCLID, 
will be able to apply these techniques to map clusters in 
much larger volumes. Such large catalogs will also be 
able to tightly constrain cosmological models (Takada & 
Bridle 2007; Dietrich & Hartlap 2010). 
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Table 5 Catalog of the Matching Groups/Clusters between v > 3.5 Convergence Peaks and 
K2-detected Groups/Clusters in CFHTLS-Wide Wl Fields. K2 redshift denotes the median 
redshift of bright (i < 20) cluster members 
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Table 6 Catalog of the Matching Groups/Clusters between i> > 3.5 Convergence Peaks and 
X-Ray/K2-detected Groups/Clusters in CFHTLS- Wide Wl Fields. K2 redshift denotes the 
median redshift of bright (i < 20) cluster members. X-ray redshift denotes the photometric 

rcdshifts of X-ray clusters (Adami et al. 2011). The parameters <J xr ay and arG tne offset 

between X-ray/optical and weak lensing center, respectively. 
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